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Abstract 

1 /3 

It is already known that a positive Cosmological Constant A sets the scale vq = (I^s^a) > 
which depending on the mass of the source, can be of astrophysical order of magnitude. This 
scale was interpreted before as the maximum distance in order to get bound orbits. In this paper 
I compute ro with a different method and obtain its first order correction due to the angular 
momentum L of the test particle moving around the source. I then re derive by using more rigorous 
methods the maximum angular momentum in order to get bound orbits L max = j^r 2 ?^) 1 / 3 and its 
corresponding saddle point position given by r x = ^(Sr^r^) 1 / 3 . Here r s = 2GM is the Schwarzschild 
radius, r\ = is the Cosmological Constant scale and /3 is a dimensionless parameter given by 

Q — L max 

P — L ■ 



PACS numbers: 



I. INTRODUCTION 



The Schwarzschild de-Sitter (S-dS) space in static coordinates has been widely studied 
in the past. Its analytic extension has been performed by Baianski and Ferrari pQ. They 
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interpreted the scale r = (Iva) as the distance where the 0-0 component of the S-dS 
metric takes a minimum value. As a consequence of this, it was found in [2J that r represents 
a transition distance after which a photon suffers a gravitational blue shift when it moves 
away from a source. The same scale is used by Bousso and Hawking in order to find the 
appropriate expression for the temperature of a black hole immersed inside a de-Sitter space 
[3]. In such a case, the distance tq is interpreted as the position of the static observer in order 
to find the appropriate normalization for the time-like Killing vector. As a consequence of 
this, there exist a minimum temperature for the black hole given by T = [31 H] . In [5] , 
Balaguera et al, found that vq represents the maximum distance within which we can find 
bound orbits solutions for a test particle moving around a source. In the same manuscript, 
the velocity bounds for a test particle inside the S-dS space were obtained, this work was 
then extended by Arraut et al in [U] in order to incorporate other metric solutions. In [SJ, 
the authors also found that there exist a maximum angular momentum L max for the test 
particle to be inside a bound orbit. If L = L max , then there exist a saddle point for the 
effective potential at the distance r x . In this paper I improve the computations given in [5] 
and then I find the first order correction for the distance ro- This correction is due to the 
angular momentum L of the test particle moving around the source. Additionally I prove 
that the saddle point position r x and as a consequence the maximum angular momentum 
L max mathematically represents an extremal condition for the discriminant of a 4th order 
polynomial. In fact, the condition L = L max is just equivalent to say that the discriminant 
for the fourth order polynomial is zero D = 0. The paper is organized as follows: In section 
[TTJ I make a review of the Schwarzschild de-Sitter metric in static coordinates, I then obtain 
the effective potential (U e ff) and analyze the validity of the coordinate system used for 



the present analysis. In section |III| I obtain the critical points of the effective potential 
^ dU ^ f = oY then I obtain the scale tq with its first order correction due to the angular 
momentum L in terms of a dimensionless parameter (3 = L "l ax . Additionally, I re derive the 



scales L max and r x . In section IV I compute by using the same methods of section III the 



zero effective potential conditions U e ff(r) = 0. The objective is to verify the consistency of 



the critical point results dU ? ff = 0. And finally in section 



V 



I conclude. 



II. THE SCHWARZSCHILD DE-SITTER METRIC IN STATIC COORDINATES 

Schwarzschild-de Sitter metric in static coordinates, is given by: 

ds 2 = -e^ r) dt 2 + e~ y{r) dr 2 + r 2 d8 2 + r 2 sin 2 Bd$ 2 (1) 

where: 
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e"M = l- r JL-* (2) 
r 3r^ 



if we want to find equations of motion for the S-dS metric, the work is simplified if we use 
its symmetry. From symmetry arguments, we know that there are 4 killing vectors, three 
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for spatial symmetry and one for time translations. Each of these will lead to a constant of 
the motion for a free particle [TJ. If is a Killing vector, we know that: 

dx^ 1 

Kfi—jr- = Constant (3) 
dX 

is a constant of motion. In addition, there is another constant of motion, namely, the 
geodesic equation (together with metric compatibility). Then the quantity: 

^ -dx -dx = ~ e (4) 



is a constant of motion along the path [7] . In eq. [4] e = 1 for massive particles and e = 
for massless particles. The quantity associated with invariance under spatial rotations is 
the angular momentum. We can think about the angular momentum as a three-vector with 
magnitude (one component) and direction (two components). Conservation of the direction 
of angular momentum means that the particle will move over a plane, we can choose this 
to be the equatorial plane of our coordinate system. If the particle is not initially over the 
plane, we can rotate the coordinate system until that condition is satisfied. Then we can 
choose the angle (TJ: 

* = \ (5) 

we then have two remaining Killing vectors corresponding to the conserved quantity 
related to time translations and the magnitude of angular momentum. The quantity related 
with time arises from the time-like Killing vector: 

K» = (d t y= (1,0,0,0) (6) 
the Killing vector whose conserved quantity is the magnitude of the angular momentum 



is: 



^ = W = (0,0,0,1) (7) 
in both cases, it is convenient to lower the index and then obtain: 



thus: 



K» = {-e^\ 0,0,0) (8) 



= (-(1-7-^) 0.0.0) (9) 



and: 



R, = (0,0,0, r 2 ) (10) 
from eq. [3j the two conserved quantities are: 



E-K^l (U) 



and: 
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r 



" d\ dr 

where L is the magnitude of the angular momentum, and r 
Developing explicitly eq. [4] for massive test particles, we obtain: 



1 = #00 



dt\ 2 ( dx l N 



dr 



+ 9 *{dr 



if we introduce the S-dS metric given by eq. [TJ then: 



3r 2 A 



dt 
dr 



this equation is just equivalent to: 



r 



3r 2 A 



(12) 

A is the proper time. 

(13) 



j) + ''\ ■ 
dr J \ get 



(14) 



r 



3r\ 



- 1 



3r 2 A 



using eqns. [2 11 and 12 we obtain: 



2 / dt\ 2 (dr 
\dr J \dr 



+ r 2 1 



3r\ 



dr 



(15) 



- 1 



3r\ 



—E 



( dr 



1--- 

r 



3r\ 



dividing by two, developing and rearranging terms, we get: 



dr 



+ 



L 2 
2? 



2r 3 



2r 



6 r\ 



E 2 + 



3^ 



C 



(16) 



(17) 



where C is a constant depending on the initial conditions of motion, we can define the 
effective potential as: 



U, 



eff\ 



2r 



then eq. 17 is equivalent to: 



lr 2 



6 r\ ^ 2r 2 



2r 3 



(18) 



dr 
dr 



+ U, 



eff\ 



E 2 + 



L 2 ^ 
3rl 



1 



C 



(19) 



the first term on the right-hand side of eq. 18 is the Newtonian gravitational potential, 
the second term is the A contribution (it reproduces a repulsive effect), the third term is the 
centrifugal force and it takes the same form in Newtonian gravity and General Relativity. 
The last term is the General Relativity correction to the effective potential. This last 
contribution is important at short scales i.e comparable to the gravitational radius 2GM. In 
this manuscript, we suume that r\ » r s . We can use a relation between L and E, obtained 
as follows: 



dr 



d<p dt 
dt dr 



r 2 d>e 



-u(r) £ 



(20) 
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where 



is the angular velocity of the test particle moving around the source. 



Replacing eq. [20] in eq. [17J we obtain: 

W I ! 



r 



4i2_-2j/(r) 



3r 2 



c 



1 

2 



(21) 



then the condition C > — | must be satisfied in agreement with [5] where the same 



analysis was performed. Additionally, eq. [17] can be written as: 

2 r2 ^,r2 „, i „J2 i/ r2 



1 /dr cftV L 2 

2 Vrftrf7y + 272 



r.L 2 
2r 3 



If. 
2r 



1 H 
6 r 2 



i 



c 



using again eq. [TT] and eq. [20] after a rearrangement of terms, we obtain: 

E 2 -- 



e 2u{r) _ 2K + </>V (1 - e 



(22) 



(23) 



r 2 + r 2 2 as twice the Kinetic energy value. This equation 



where we have defined 2K 
offers certain restrictions, since we need E 2 to be positive definite. The condition E=0 
is equivalent to the condition e u ^ = 0, namely, the event horizon condition. Then the 



expression 23 shows the regime of validity of our coordinate system. The event horizons in 
S-dS metric, are obtained under the condition: 







1_1±_ 
2r 



3r 2 



it is already known that the solution of eq. 24 provides two real solutions given by: 



>'< 7/ -2r A ro.s [ i | cos 



and: 



r BH = -2r A cos 



cos 



-i 



3r s 
2r A 

3r s 
2r A 



2tt 



(24) 



(25) 



47T 



where tch is the Cosmo logical Horizon and tbh is the Black Hole event horizon. In this 
paper I will work in the region between the two horizons. I also assume r\ » r s . If r\ ~ r s , 
then Tch ~ t bh and then our coordinate system is inappropriate [3]. I will not consider 
such situations in this manuscript. Under the condition » r s , eqns. 25 become [2]: 



r C H 



\/3r A 



(26) 



and: 



tbh ~ r s + 



lrj 
6 r\ 



for completeness, the angular momentum in terms of the variables introduced in eq. 23 



is: 
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r 402 e 3„(r) 

' e 2u(r) _ 2K + 0V 2 (1 - e^ r ); 
as can be easily verified. Additionally, the constant C, can be expressed as: 



(27) 



C(ri,Vi,<j)i) 



v? + He Mn) +v}i)Ueff{ri)L=0 + ^< 



Sri 



2 (e 2!/ ( ri ) 



(28) 



where U, 



eff{ri) L=0 



r s 

2n 



l t 
6rf> 



and the subindex z denotes initial conditions, v is the 
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as L — > 0, 



total velocity of a test particle, and = rcf) is its tangential velocity. Eq. 
recovers the result obtained in ref. [5], where the same analysis for the case L = was 
performed. 



III. CIRCULAR ORBIT CONDITIONS FOR THE EFFECTIVE POTENTIAL 
WITH A 



We can know obtain the critical points for the effective potential obtained from eq. [T8j I 
will start the analysis with the case where the Cosmological Constant goes to zero (A — 0). 
We know that the effective potential for a massive test particle is given by [7]: 

r, L 2 r q L 2 , . 

W"--^^— (29) 

taking the derivative with respect to r of this equation we then obtain: 

dU e ff(r) A=0 _ _ L? 3r s L 2 

dr ~ 2r 2 r 3 2r 4 1 ' 



in order to find the circular orbits, we have to send to zero eq. 30, then: 



r 2 - —r c + 6L 2 = (31) 

where we have defined r c as the distance at which the orbits for a test particle moving 
around the source are circular. The solutions of this previous equation are: 



L 2 ± JL* - 3r 2 L 2 , . 
^ — (32) 



if L — > oo, then: 

L 2 ±L 2 (l- 3 S) f2L 2 3 



• 2 r. ) (33) 



the first root is the stable orbit [7] and the second root is proportional to the Schwarzschild 
radius and it is the unstable equilibrium orbit. The minimum angular momentum in order 



to get bound orbits is obtained as the discriminant of eq. 32 is zero, in such a case 



Lmin = v 7 ^ (34) 
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when the angular momentum takes this value, the two distances for circular orbits ob- 



tained in eq. 32 become the same and they correspond to the saddle point position given 
by: 



r cx = 3r s (35) 

the same result is obtained if the first and second derivative of the effective potential 
given by eq. [29] go to zero (See appendix [A]) . If we replace the second root of given in eq. 



33 



,) inside eq. 



30 



then we obtain: 



dU, 



eff\ 



dr 



A=0 



2 

9^ 



8L 2 8L 2 
+ 



27r 3 27r 3 



2 

9^ 



(36) 

This distance becomes 
such a case, then the 

must be ignored from eq. |30 at short distances r. By the same argument, the term 

Then for short scales 



we then conclude that r r 



is not an exact solution for eq 



to be an exact solution only under the approximation L » r 
term L 

— — must be dropped from eq. 



• x 7ln 



A = can be approximated: 



29 



r — > ?% 



U ef f{r) A=0 



17_ 
2^2 



r s L 2 
2r 3 



the effective potential as 



(37) 



ignoring then the standard Newtonian contribution. As a proof of this statement, we can 
take the derivative with respect to r from the previous equation: 



dUeffjr) 
dr a=o 

if we set this result to zero, then we obtain: 



L 2 3r s L 2 

I — 

r 3 2r 4 



(38) 



(39) 



which is just the second root given in eq. 33 Fig. [T] shows the effective potential curve 
for distance scales around the gravitational radius GM. If we accept the fact that the term 
— is negligible for short distances r in the equation [29] as far as the approximation L » r s 
is valid, then the same argument will be valid as 0. And the additional contribution 
to the effective potential given by — will be also negligible for short distances r. Later 
we will see that the A effects become important for distances of the order of magnitude of 

r ° = (§ r s r A) 01 l ar g er - The full effective potential energy with the inclusion of the A 
term is given by: 



U, 



eff 



T = — 



2r 



Qr 2 , 2r 2 



r s L 2 
2r 3 



(40) 



under the previous approximations, the first circular orbit which is the unstable one will 
be located at the same distance r\ = |r s , even for the case A / as far as the condition 
rA >> L » r s is satisfied. Figure [T] shows the effective potential for distance scales of 
the order of magnitude of the gravitational radius r s . If we consider the second position for 
circular orbits given by r 2 = — (see eq. 
not an exact solution of eq. [30] neither. In 
we then obtain: 



33|) as A = 0, it is then possible to verify that it is 
act, by replacing r 2 in the mentioned equation, 
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dr 



A=0 



3r| 
32L 6 



(41) 



then, r 2 is not an exact root for the case A = 0. It is easy to demonstrate that r 2 becomes 

is neglected. Then 



30 



an exact solution only if the term 3 t^4 (GR contribution) from eq. 
it is necessary to ignore the term — ^fpr in eq. 29 in such a case. Then we conclude that for 
large distances r, as A = 0, under the assumption L » r s , the effective potential (A = 0) 
can be approximated as: 



Ueff(r)A=o~-^ + ^r 2 (42) 
if we take the derivative of this expression with respect to r and send it to zero, we get: 



dU, 



eff\ 



2r 2 







^ , , (43) 
dr a=o 2r z r 6 

then the solution for this equation is of course r 2 as it should be. This result will be 



important at the moment of making the computations for the case A ^ 0. From eq. 40 



we see that the A term should be important for large values of r. If we still satisfy the 
condition L » r s , then the effective potential for the case A ^ and for large scales can 
be approximated like: 



U, 



eff 



r 



2r 



Qr 2 A 



I?_ 
2^2 



(44) 



where we have neglected the general Relativity correction term in agreement with 

the previous assumptions. If we take the derivative with respect to r of the previous equation, 
we obtain: 



dU, 



eff\ 







dr 2r 2 3r A r 3 

multiplying this equation with r 3 and rearranging terms, we get: 



(45) 



2 T 2 



r — -r A r s r + 3r A L 







(46) 



this equation is a reduced fourth order polynomial in agreement with appendix [Bj We 
can make a comparison with the standard reduced form given by equation B7 of appendix 
[B| Here I rewrite the equation in terms of the variable r as: 



r 4 + Pr 2 + Qr + K = 



if we compare this equation with the result [46j then we can define the coefficients: 

3 o 



P = Q 



-r A r s 



K = 3riL 



2 r 2 



(47) 



(48) 



we can extract the associated cubic equation. This is given by eq. B9 , which in this case 
becomes: 



z 3 - \2r 2 A L 2 z - -Art = 



(49) 
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in agreement with the coefficients defined in eq. 48 The solutions for the fourth-order 
equation 46 in terms of the associated solutions of eq. 49, can be found from eq. BIO which 
here is just rewritten in terms of the variable r as: 













* 

r x = 


I, 








r* = 


I, 


fz~\ - ^ 






* 1 










3 2 










* 1 










r4 = 2 


(-> 


fz[- . 


/Z2 ~ 1 





(50) 

also the solutions of the associated third-order equation must satisfy the condition |B11| 
from appendix [Bj Q has been defined in eq. 47 and in this particular case is given by the 
result 48 Our objective here is to find the real roots for the fourth-order equation. We can 



49 which is already in a reduced form, with the standardized reduced form 



compare eq. 

given in eq. C5] from appendix O which is rewritten here in terms of the variable z as: 



z + pz + q = 



where p, and q are constant coefficients. In our case, from eq. 49 we obtain: 
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V 



■\2r\L 2 



q 



4 2 

~J A r s 



(51) 



(52) 



it is necessary to find the discriminant D which is given in eq. C7 rewritten here as: 

D = ^P 3 + -/ (53) 



where p and q are given by eq. 52, replacing, we have: 

51 



D 



-64r 6 A L 6 + 



8^4 



64 



-r A r 



A' s 



(54) 



the solution for the associated cubic equation given in 51 , depends on the sign of the 



discriminant defined in eq. 54 Can be demonstrated that if D < 0, then there are no real 



solutions for the fourth-order equation 46 This case then is not important for us. However 
if D > 0, we can get real solutions. Then D = must provide some limit condition given 
by: 







-64r b A L b + —r s A r 



if we solve for L, we obtain: 



3 2/3 



64 



{r 2 s r A ) 



1/3 



(55) 



(56) 



this condition is coincident with a result obtained in [5]. In such a case, it was interpreted 
as the maximum angular momentum for a test particle to be inside a bound orbit. Here that 
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interpretation is verified by the fact that for the case L > L max there are no real solutions 
for the fourth order polynomial. We can now concentrate on the case D > 0. We also know 



from eq. 52 that the conditions p < and q < are satisfied. If we want to find the solution 



for the associated third-order equation 51, then we need to define first the parameter given 



in eq. C8 of appendix [Cj In this particular case that parameter R depends on the angular 
momentum and it is given by: 



R = -2r A L 



(57) 



here we use the case ii) of appendix [Cj specifically by eq. Cll From eq. 52 and eq. 57 
we obtain the auxiliary angle: 



solving for <p gives: 



coshcj) 



cosh 



9 tatI 
64 T 3 " 



l ( 9 r A r\ 



64 L 3 



if we introduce eq. 56 in this result, then the auxiliary angle can be expressed as: 



cosh 



i ( L 



L 



(58) 



(59) 



(60) 



the three solutions for the associated cubic equation 49, are obtained from eq. C12 of 
appendix O However in this case y = z, so we get for the first root: 



Z\ = —2Rcosh 



replacing the results 57 and 60, we get: 



(61) 



:.i = !/■ \Lcosli | -cosh 1 (~^~ 



the second root of the associated cubic equation 49 



is: 



z 2 = Rcosh (J^j + iVSRsinh (J^J 



replacing again the results 57 and 60 we obtain: 



(62) 



(63) 



Z2 



-2rALcosh ^-cosh 1 ^ max ^j \ _ 2%\[?>r ALsinh ^-cosh 1 



finally, the third root is given by: 



^3 = z* 2 = Rcosh 
thus by analogous procedure, we get: 



i\/3Rsinh 



(64) 



(65) 
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: :; = :.; = — 2r \Lcosli ( -cosh 1 ( ™ x 



+ 2iy/3r /iLsinh f -cosh 



1 / -L-'max 



(66) 



we can define the dimensionless parameter L = L ™ ax . Then eqns. 
equivalent to: 



62 



64 



and 



Zl = 4 



P V o 



max 

z> = -2—p-cos.. Vi 



h ( lcosh-\pf) - 2iV3 T -^P^sinh ( ^-cosh- l {Pf 



and: 



Z3 



p \3 / 



^^sinh I -cosh 1 ((3) 3 
P V o 



if /3 = 1, then L = L max and: 



zi = Ar A L r 



66 



arc 



(67) 
(68) 

(69) 
(70) 



z 2 = -2r A L r 



*3 



(71) 



additionally, the solutions must satisfy the condition |B11[ If we take the root square of 
eqns. 



70 and 71 Then the condition Bll becomes 



2 \/ T A^max \j 2t ' A^max \J 2f ' A^max — Q 



-r*r, 



(72) 



the left hand side of this equation becomes \/ —2r j\L max = iy/2rj\L max . Then eq. 72 

3 „ 



is: 



- 4(r yl L, 



^3/2 



s> A 



(73) 



which is consistent with the maximum value of the angular momentum found in eq. [56 



If we replace the results 70 and 71 in eq. 50, then there is only one real solution for the 
fourth order equation given by: 



r x 



r l = r 2 = V r A^max = T x (74) 

where r x is the the saddle point position for astrophysical scales. If we use eq. [56] in the 
previous solutions, then we get: 



r, = r = r x 



(3) 



1/3 



-(r s r A ) 



2x1/3 



(75) 



which is the same result obtained in [3] by using different methods. We will now expand 
cosh~ l ((3 3 ) by using the approximation: 



cosh 1 x = In (^x + V x 2 — 1 



(76) 
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if x 2 » 1, then: 

cosh~ l x w ln(2x) (77) 
in our particular = (5 3 . Then we can use the approximation as /3 6 >> 1: 

cosh-\(3 3 ) = ln{2(3 3 ) (78) 

then we can study easily the cases for which the condition /3 6 >> 1 is satisfied. This 
condition is equivalent to say that L 3 max » L 3 . Then we can approximate eq. 62 like: 



thus: 



2i ~ /^f^cosh {Ln{2fl 3 p) (79) 



z 1 ^2r A L max ^ 3 + ^^j (80) 



the root square of this result is: 



v^=^/2r^ ma ,((2)V3 + _L_^ (81) 
under the condition (3 2 » 1, an expansion around ^ can be performed and we get: 

: j = 4 1/3 a/ r A L max + ^— a/ r A L max (82) 



doing the same for z 2 given by eq. 64 we obtain: 



Z2 ~ _ 2 !^^ C0S / l (Ln(2) 1 / 3 /3 - 2y/Zi r ^^sinh{Ln{{2) 1 l*P) (83) 

then: 

z, * -r A L max ((2)V3 + _L_) _ y^L^ ((2)^3 _ _L_) (84 ) 
if we take the root square for this result, then: 

- ^r A L max ((2)V3 + _L_) + V3 irA L max ((2)V3 _ _L_) (85 ) 
this result can be translated into a complez number in polar representation given by: 



yfr ~ M /2(2)V3 r ^_ e W3 + ^ ^ e-Ws (86) 



(2)2/3r A L 

Expanding again in a Taylor series, we get: 



% * i U*yfcj^»e*l* + l^ e -/ 2 ) (87 ) 



2 /3 2 
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by using the same procedure for z 3 : 



4: 1/3 ^/r A L max e~ 



+ 



lvCT 



/? 2 



(88) 



where we used the condition z 3 = z\. Then, using again the result 50, we get as a real 
solutions: 



3 1 

IJ 2 



(89) 



inserting the result 56 in this previous equation, we obtain: 

2L 2 



(90) 
's 

this result is just the standard circular stable orbit obtained even for the case A = 0. 
Fig. [2] illustrate the behavior of the effective potential for intermediate distances scales. The 
stable orbit condition is clearly perceived from the curve. Here we are working under the 
condition /3 2 >> 1 is satisfied. Additionally, there is a second real solution for the equation 



46 This is r* 2 in eq. 50 given explicitly as: 



4 1/3 \/r A L max - — a/ r A L r . 



using again eq. 



56 , we obtain: 



-rjr 



s< A 



1/3 



;(3r s r 



2 W3 



(91) 



(92) 



this distance can be of astrophysical order of magnitude. Fig. [3] illustrates the behavior 
of the effective potential for distances near r' . In this figure it is clear the departure with 
respect to the behavior found in the standard Schwarzschild case. If 1 » 1, then: 



r 2 w 



-r,r 



s' a 



1/3 



(93) 



consistent with the results obtained in [5] as L — > 0. On the other hand, it is possible to 
obtain an approximate relation between r\ and r|. If we replace eq. 56 in eq. 89, then we 
get: 



if we solve for (3, then: 



i3 2 



if we replace this in eq. 92 , then: 



(3r s r2)V3 



:r K r, 



1/3 



-r-i 
3 1 



(94) 



(95) 



(96) 
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this is an approximate relation among r* and r\ which correspond to circular orbits, one 
stable (r*) and the other unstable r\ and due to a local effect of A. This previous relation is 
valid under the condition L » r s . Figure [2] shows the position for the stable circular orbit 
and E] shows the same but for the unstable one. 



IV. DISTANCES FOR ZERO ENERGY EFFECTIVE POTENTIAL 



I will study the scales for which the effective potential given in eq. 40 is just zero. The 
procedure will be analogous to that performed in the previous section, taking the convenient 
approximation when we are at short distances or at large distances from the source. The 
general equation to be solved now is: 



U ef f(r) = - 



LP 



r,L 2 



2r Qr\ 2r 2 2r 3 







(97) 



again for short distances r relative to the source, if L » r s then eq. 37 is valid and the 
previous condition reduces to: 



Ueff(r) 



rX 2 



2r 2 2r 3 

this equation is easy to solve, the result is simply: 







(98) 



•7 

Ti = 7%. 



i - ,. (99) 

where the super index means Zero. Figure [T] shows the plots of the effective potential for 
short distances in comparison with r . Different plots correspond to different values for the 
angular momentum. 



On the other hand, for large distances relative to the source, we again use eq. [44] and 
then obtain: 



^ (r) "-27-6^ + 2^ = ° 



(100) 



this equation can be transformed as: 



r 4 + 3r\r s r — 3r^L 2 = 



(101) 



working with the same procedure used in the last section, from equation [47] of section [ITT 
we get: 



P = ; Q = 3r 2 A r s 
from eq. B9 in appendix [Bj we have: 



K = -3r 2 A L 2 



z 3 + \2r\L 2 z - 9r A A r 2 = 



comparing this equation with eq. 51, we have: 



2 r 2 



p = \2r\L 



-9r> 2 



(102) 



(103) 



(104) 



then from equation 53, we have: 
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FIG. 1: Effective potential for distances of the order of the gravitational radius r s and L >> r s 
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D = 6Ar b A L b + —r\rl > 



(105) 



from eq. C8 , we obtain: 



R=- x Ar\L 2 = —2r\L 



(106) 

for D > and p > which is the case here. From appendix [Cj we need to define the 
auxiliary angle to be: 



sinhcf) = 



2R 3 



if we replace eq. 104 and eq. 106 in this previous equation, then: 



sinhcf) 



9 tatI 



from eq. 56 , we already know: 



(107) 



(108) 



it is convenient to write this equation as: 



(109) 



4L: 
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replacing this in eq. 108 



(110) 



sinhcf) = 4 



if we use the definition (3 = then: 



sinhcf) = A/3 3 



and then solving for /3, we get: 



(111) 



(112) 



= sink' 1 (40*) 



(113) 



using now equations |C14| with z = y, and introducing eq. |106| and eq. |113| we obtain the 
solutions for the associated cubic equation: 



Z\ = Ata m ^ x sinh ( -sinh 1 (4/3 3 ) 
P V'J 



(114) 



z 2 = —2ta ™l ax sinh ( -sinh 1 (4/3 3 ) 
P \3 



2V3ir A I ^cosh (\smh~\Ap\ 
P \<j 



;ii5) 



and: 



Z3 



-2r A ^sinh Q S m/i- 1 (4/3 3 )^) + 2VSir A ^cosh Qszn/T 1 ^/? 3 )^ (116) 
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following the procedures of section III I will analyze two cases: 



Case i). (3 = 1 

For this case, equation |1 14| becomes: 



thus: 



Zi=±r A L max sinh[-sinh-\4) 



the same procedure on eq. 115 and we obtain: 



22 



which is approximately: 



or in polar coordinates: 



2 2 



2 2 ~ ~2 rALr 



21e i(? ' /5)n r A L 



- ,V3ir A L 



21e- i{6/w)n r A L 



eq. 



116 then becomes: 



3 T 5 r T 

23 ~ --r A L max + -V3ir A L r , 



or in polar coordinates: 



2.3 



21 e -*(7/5) V r 



2i e *(6/io) v r 



;ii7) 



(118) 



-2r A L max sinh Qsm/i x (4)^j - 2v / 3«r /1 L maa .cos/i Qsm/i x (4)^j (119) 



(120) 



121) 



(122) 



(123) 



the three solutions must satisfy the condition Bll In fact, by replacing eqns. 118 121 



and 123 inside the mentioned condition, we get: 



^3r A L max 2{2l llA )^Jr A L max ^Jr A L max « 3r 



s> A 



(124) 



after replacement of eq. 56, this result is in agreement with eq. 102 Then the roots 



corresponding to the fourth-order equation 101 are given in this case by eqns. 50 , but only 
the last of these group of equations gives us a real and non-negative solution as can be easily 



checked. If we replace eqns. 118, 121 and 123 inside eq. 50, then 



I (-V^a~L^ x - (21)^^J^ x e^^ - (21)VV^ 10 >) (125) 



by using the identity e = cosO + isin9, then 



^ ((-\j3r A L max - 2(21) 1/4 ^/r A L max cos (J^j 



(126) 



or: 



0,4:^/r A L max ~ 4x10 1 y / r A L r 



(127) 
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from eqns 74 and 75 , the previous equation can be written like: 



4x10" 



(128) 



where r x symbolizes the saddle point position previously defined. Then the effective 
potential goes to zero before the saddle point position as it must be. Then at the saddle 
point position, the effective potential is negative. In order to verify it, we can evaluate the 



effective potential 40 at the saddle point position 75 and we then obtain: 



U, 



r s 
2r„ 



lr 2 
— - + 
6 r\ 



2r2 



2rl 



(129) 



but at the saddle point, the distances are large enough and we can neglect the GR 
correction term as far as the condition L » r s is satisfied. Then if we consider eq. 44 



after replacing the result 75 then 



U, 



eff\ 



2r„ 



6r 2 A 2r 2 



2/3 



< 



(130) 



A *'x \ r Aj 

then the effective potential at the saddle point position is negative under the previous 
assumptions. Although it is in general very small in magnitude since the term ^ << 1, in 
fact this term is of order 10~ 15 asM~ M sun (the solar mass), however its value depends of 



the mass of the source producing the gravitational attraction. Eq. |128| can be expressed in 
terms of the maximum angular momentum given in eq. 



56, then: 



r, ~ 



r2 

max 



(131) 



Fig. |2j shows the zero for the effective potential located at scales of the order of magnitude 
of r . From Fig. |3j it is clear that there are no zeros for distances larger then r . For f3 = 1 
we know at what positions the effective potential goes to zero (see eqns. 99 and 128), we 
know the position for the saddle point (see eq. 75). We know at what position the effective 
potential takes a maximum value (see eq. 39) and we we also know what is the value for 
the effective potential at the saddle point (see eq. 130). For a final picture of the important 



points for the effective potential curve, we need to find the value for the effective potential 
when it is evaluated at the distance given by eq. [39l This distance corresponds to a local 



maximum and is located at short distances r with respect to the source. Then we need to 



use the approximation given by eq. 37, which evaluated at the distance 39 gives 



U, 



eff 



2 / L 
2' 7 ^~27V^ 
at the maximum angular momentum value |56l we get: 



V, 



eff \ 



2/3 



216 \r s J 



> 



(132) 



(133) 



which is in general much bigger than one. 



Case ii). (3 > 1 



18 



c effect 




c effect 




t'effeet 



^effect 




f effect 



-1.x 10 


-5 




-2.x 10 


-9 




-3.x 10 


-9 




-4.x 10 


-9 




-5.x 10 


-9 




-6.x 10 


-9 






^effect 





2.x 10 7 4.x 10 7 6.x 10 7 8.x 10 7 1.x 10 s 




1 erect 




FIG. 2: Effective potential for scales r s « r « ro and L >> r s . 



Let's study now the general case for any (5 » 1. Then we can use the approximation 
sinh~ 1 (4/3 3 ) rs Ln(8/3 3 ). With this in mind, inserting this approximation in eqns. 114, 115 



and 116, we obtain: 



Z\ ^ 4 



n 



'-sink (Ln(8/3 3 ) 1/3 ) 



(134) 



thus: 
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FIG. 3: Effective potential scales of the order r ~ r$ and L » r s 



zi « 4r /1 L maz ( 1 - — ) (135) 



4/3 



by taking the root square, then we get: 



z x « 2\lr A L max ( 1 - ^ ) (136) 
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for (3 2 » 1, we can expand in a series: 



for 2 2 and z 3 we can use analogous procedures: 



1 

ma. r 1 1 TTTlO 

8p 2 



(137) 



z 2 « -2 







sira/i (Ln^ 3 ) 1 / 3 ) - 2x/3z^p COS /i (Ln^ 3 ) 1 / 3 ) (138) 



/3 



or: 



z 2 ~ -2r A L 



1 

1 ~~ i/3 2 



for ^3, we have: 



:» * -2r,Z,„„, ( 1 ^ 



2\/3iryiX mox ^1 + -^J 
2v / 3zr yl L m(M . ^1 + 



(139) 



(140) 



using the same reasoning done in order to obtain eq. 124 for the case i), as the cases 
must have compatibility with each other, then we can express directly the complex numbers 
given in eq. |139| and eq. |140| as: 



z 2 « 2r A L max (-1 - + 2^^ (l - V3i) 



or: 



and: 



4/3 2 



7n ~ Ar a T P - i27T / 3 4- rALmax -fa/3 

Z 2 ~ yl-^maxC + 



/3 2 



(141) 



(142) 



z 3 « 4r /1 L_e^/ 3 + I^^ 3 



taking the root square of eq. 142 and eq. 143 we get 



(143) 



l Ar A L max e-^lz + 



r A L r 



& 1 



and: 



l 4r A L max e^ + r -^^ ( ' 



if we execute the expansion in a Taylor series, for large values of /3 2 , we get: 



and: 



(144) 



(145) 



(146) 



z 3 « 2A/r j4 L ma:r e J7r/3 + — y^A^r 



(147) 
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replacing eqns. 137 146 and 147 inside of 50, the only physically acceptable real solution 
is given by: 



z 3 \fr~A~L 

max 

r> — 

3 8 P 2 

introducing eq. 56 inside of this result and using the definition L = we get: 



148) 



(149) 



if we compare this result with the one obtained in eq. [90J it is clear that the effective 
potential goes to zero before the second critical value (circular orbit condition) which is 



completely consistent. Then at the value given by eq. [90j the effective potential must be 
negative. There is no other real zero for the effective potential in this second case and at 



big values of r, which means that at the third critical point given by eq. 92 (which is a 



maximum), the effective potential is still negative. If we evaluate the effective potential eq. 



40 at the distance 90, term by term, we get: 



V, 



off 



21/ 



+ '- 



1 (1 



2/3 



2/3 



2/3 



/3 2 - 



i a 



6 V4 



4/3 



V2rJ /3 4 



2/3 



^ v 3 y \r A ) 

for most of the situations (3 » 1, then: 



/3 2 - 16 



8/3 



rA 



4/3 



& 1 



(150) 



U, 



eff 



2L' 



1 (1 



2/3 



r_s\ 
TAJ 



2/3 



/3 2 - 16 



8/3 



4/3 



/3 4 <0 



(151) 



r s J 2 \3J \r A J '' " \3 J \r A J 
the last term in this equation can be neglected. If we want the two terms on the right-hand 



side of eq. 151 to be comparable it is necessary: 



i n\ 2/3 /r i y /s 



2 V3 



8/3 



4/3 



solving for (3, then: 



I? 



(32) (rj 



2/3 



(152) 



(153) 



for a source with a mass of the order of our sun: 



w 10 7 



(154) 



which is outside of the possible range of values taken by (3 in order to keep the assumption 



L » r s as can be easily verified. Then equation 151 can be approximated as 



U, 



eff 



2L' 



1 /2 



2 V3 



2/3 



ta) 



2/3 



(3 2 <0 



(155) 
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if we evaluate now the effective potential at the third critical point given in eq. 92, then 
we obtain term by term using the approximation 44 



off 



1/3" 



3V3 V 2r A 



2 /3 gl/3 



6 \2r A 



2/3 



+ 



1 /3\ 4/3 (rgryQ 2 / 3 
(6r s r A ) 2 / 3 



2/3 2 



(156) 



where we have used the definition L = and 



mately: 



56 



for big values of j3, we get approxi- 



U, 



((Sr s r A ) 1/3 ) - " (^) 



2/3 



< 



(157) 



which is equal in order of magnitude to U e ff evaluated at the saddle point given in eq. 



75 (see also eq. 130). 



V. CONCLUSIONS 

In this manuscript I could re derive some of the astrophysical scales (effects) due to the 

I/O 

presence of a Cosmological Constant A. The scale r = (§r s r^J was derived again by 
using a more rigorous method in order to get the first order corrections due to the angular 
momentum of a test particle moving around a source. By taking into account its first order 

1 /3 

correction, it takes the value r* 2 = r' ~ (| r s r i) ~ (3r s r^) 1//3 , where (3 = ^^a. This 
is the new contribution of the manuscript. Additionally, I demonstrated that the condition 
for getting the saddle point position r x is just equivalent to setting to zero the Discriminant 
D = for the fourth order polynomial. From this limit condition, the maximum angular 
momentum L max in order to get bound orbits is obtained. There is a saddle point position 
r x for the effective potential when the test particle has an angular momentum L max . 
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Appendix A: Minimum angular momentum with A = 



In the effective potential given by eq. 29 , we can find the first derivative, which is given 



by eq. [30j also, we can find from this the second derivative given by: 



d 2 U ptt (r) 



dr 2 



r s L 2 r«L 2 
+ 3— -6-V- = 



(Al) 



where we have sent to zero the second derivative of the effective potential; this equation 
can be rewritten as: 



r 2 - 3 



r + 6L 2 = 



(A2) 



the equality between this equation and that gotten in 31 , for which we got in that moment 
the condition for circular orbits, implies: 



or2 

2 i n T 2 2 

r r + 3L = r — 



3 (2L 2 



2 V r„ 



r + 6L 2 



solving for r, we get: 



r = 3r„ 



replacing this result inside eq. A2, we get: 



9r 2 - 9L 2 + 6L 2 = 



solving for L, we get: 



V3r 



■"mm v OT.q 

which is the minimum angular momentum. 



(A3) 

(A4) 
(A5) 
(A6) 



Appendix B: General solution of a fourth-order polynomial 

The standard form of a fourth-order equation can be taken to be: 

Ax 4 + Bx 3 + Cx 2 + Dx + E = (Bl) 
to get the reduced form of this equation, we need to do the following variable change: 

B 



solving for x, we get: 



y = x 



x = y 



AA 



B 
AA 



(B2) 



(B3) 



replacing this result inside of Bl we get: 



B Y 



B X 



B X 



Mv-u) +B [»-iA) +c \ y -fA) +d {v-Ia) +e = (B4 > 
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developing parenthesis, we get: 



A I y 4 / 5 V.3 , 3 ( B \\.2 1 ^ 



c I y 2 - 



1 (B 



2A' y 16 \ A 
B \ 



+ D[y-—j+E = (B5) 



regrouping common factors and dividing by A, we get: 



CP P 



5 



2A 2 AT 256 V A 
the standard reduced fourth-order equation is given by: 

y 4 + Py 2 + Qy + K = 



1 C 9 DB E 
16 A 3 AA 2 A 

(B6) 

(B7) 



comparing this form with B6 , we get 



K 



E 
A 



Q 
DB 



P 



1 (B 



+ 



8 \A 
1 OB 2 



4A 2 16 A 3 



c 




B 


A 


-K 


A 


3 


CB 






2A 2 


+ 




-( 






256 I 


^4 



D 

A 

d 



the associated cubic equation of |B7| is: 

z 3 + 2Pz 2 + (P 2 - AK)z - Q 2 







(B8) 



(B9) 



where P, Q and K are given in B8 Given the solutions of the associated third-order 



equation B9, the solutions for the reduced equation B7 are given by: 



2/i = g (v^i + v 7 ^ - v^a) 



Z2 + V Z 3 



^2 - V z 3 



(BIO) 
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also, the solutions of the associated third-order equation |B9[ must satisfy the following: 

Q = V^V^V^ (Bii) 



Appendix C: General solution of a third-order polynomy 

The standard form of a third-order polynomy is given by: 

ax 3 + bx 2 + cx + d = 
the normal form is gotten dividing by a this result as follows: 



where obviously, we have: 



x + rx + sx + t = 



bed 
r = - s = - t = - 
a a a 



(CI) 



(C2) 



(C3) 



with the extra condition a / 0. The reduced form of the third-order equation C2 , requires 
the change of variable: 



y = x + 



and the reduced form is given by: 



y +py + q = 
with the corresponding coefficients given by: 

r 2 2 o rs 

p = s q = — r hi 

v 3 H 27 3 



(C4) 



(C5) 



(C6) 



where r, s and t are given in C3. On the other hand, it is necessary to establish a 



classification criteria for the reduced third-order equation given in C5 , the criteria is based 
in a parameter D given by: 

also for the same classification, it is necessary to find an auxiliary parameter given by: 



R = sign(q)' 



(C8) 



where q and p are defined in C5 and C6 On the other hand, we define an auxiliary 
angle 0, which is defined depending of the following cases: 



Case i). p < 0,D < 0. 

In this case, the auxiliary angle is defined as: 
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cosct> = ^ (C9) 



with the corresponding solutions: 

Hi = —2Rcos 
3 3 

> 47 

3 3 J 

Case ii). p < 0, D > 0. 

In this case, the auxiliary angle is defined as: 



3 

y 2 = —2Rcos — + 



y 3 = -2Rcos I £ + ^ ) 



2i? 3 

and the corresponding solutions are: 



yi = —2Rcosh— 
3 

■wo = Rcosh— + iVSRsinh— 
y 3 3 

2/3 = 2/2* = Rcosh^- — iV3Rsinh^ 



Case iii). p > 0, D > 0. 

In this section, we define the auxiliary angle to be: 

Q 
2R 3 



with the corresponding solutions: 

yi = —2Rsinh 



3 



j/2 = Rsinh— + iV^Rcosh— 
o o 

2/3 = 2/2* = Rsinh— — iVSRcosh— 



(CIO) 



cosh(f)=^- (Cll) 



(C12) 



smA0 = (C13) 



(C14) 



In the mentioned three cases, D is defined in equation C7 p is defined in equations C5 



and C6 , and the parameter R is defined in equation C8 
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